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(5)  INTRODUCTION 

Of  the  genetic  alterations  associated  with  breast  cancer,  changes  in  p53  are  the  most  frequent  and 
identified  in  20-40%  of  all  cases  (Borresen-Dale,  2003;  Ziyaie  et  al.,  2000).  In  fact, 
approximately  half  of  the  major  forms  of  cancer  contain  p53  mutations,  and  the  vast  majority  of 
these  cluster  in  conserved  regions  or  “hot  spots”  (Hainaut  and  Hollstein,  2000).  Missense 
mutations  leading  to  amino  acid  changes  are  the  most  common  p53  alterations  in  breast  cancer, 
as  in  other  tumors  (Hainaut  and  Hollstein,  2000).  Together,  these  observations  suggest  a 
requirement  for  a  putative  oncogenic  contribution  conferred  by  many  TP53  mutations  in  breast 
cancer,  and  imply  that  the  development  of  small  molecule  compounds  that  may  bind  and 
reactivate  the  protein  product  of  tumor-derived  TP53  mutations  may  have  therapeutic  use  for  the 
treatment  of  breast  cancer. 

The  TP53  gene  encodes  the  p53  protein  that  regulates  the  transcription  of  a  number  of  genes 
involved  in  cell-cycle  arrest  and  induction  of  apoptosis  in  response  to  cellular  or  genotoxic  stress 
such  as  DNA  damage  or  hypoxia  (Bargonetti  and  Manfredi,  2002).  The  transcriptional  activity 
of  p53  is  mediated  by  a  tetrameric  form  of  the  protein  that  binds  DNA  in  a  sequence- specific 
fashion  to  activate  or  repress  the  transcription  of  target  genes  (El-Deiry  et  al.,  1993;  Friedman  et 
al.,  1993;  Halazonetis  and  Kandil,  1993;  Stenger  et  al.,  1994).  p53  contains  four  functionally 
distinct  domains:  a  N-terminal  transcriptional  activation  domain  (residues  1  to  44),  a  central  core 
(residues  102  to  292)  containing  a  DNA  binding  domain,  a  tetramerization  region  (residues  320 
to  356),  and  a  regulatory  domain  (residues  356-393)  (Cho  et  al.,  1994;  Pavletich  et  al.,  1993; 
Wang  et  al.,  1993).  The  vast  majority  of  tumor-derived  p53  mutations  are  localized  to  the  p53 
core  domain  (Cho  et  al.,  1994).  The  X-ray  crystal  structure  of  the  monomeric  core  domain  of 
p53  bound  to  DNA  has  provided  invaluable  insights  into  how  several  tumor-derived  mutations  in 
p53  disrupt  its  activity  (Cho  et  al.,  1994).  Specifically,  these  studies  reveal  that  the  tumor- 
derived  p53  mutations  that  are  localized  to  the  core  domain  result  in  two  different  classes  of  p53 
protein  alterations:  (1)  reduced  protein  thermostability  mutations  and  (2)  mutations  that  directly 
disrupt  protein-DNA  contacts.  Both  classes  of  mutations  functionally  compromise  the  ability  of 
p53  to  carry  out  its  normal  tumor  suppression  function  and  thus  contribute  to  neoplasia.  The 
goal  of  our  studies  is  to  identify  lead  compounds  that  bind  and  stabilize  the  subset  of  tumor- 
derived  stability  mutants  within  the  p53  core  domain.  We  anticipate  that  the  identification  of 
such  compounds  will  serve  as  a  scaffold  for  the  preparation  of  small  molecule  drugs  for  the 
treatment  of  p5  3 -mediated  breast  cancer. 

The  Specific  Aims  of  the  proposal  are  to  (1)  Determine  the  high  resolution  X-ray  crystal 
structure  of  the  p53-core  domain  bound  to  a  stabilizing  peptide  called  FL-CDB3,  (2)Use  the 
Multiple  Solvent  Crystal  Structures  (MSCS)  technique,  to  identify  novel  p53  stabilization  sites, 
(3)  Use  the  structural  information  of  aims  1  and  2  as  a  scaffold  for  using  computational  strategies 
for  the  further  development  of  small  molecule  compounds  and  peptides  for  the  reactivation  of 
tumor  derived  p53  mutants,  and  (4)  Functionally  characterize  the  p5 3 -stabilizing  and  p53- 
reactivation  properties  of  the  molecules  derived  from  aim  3,  and  determine  their  structures  in 
complex  with  p53. 
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(6)  BODY 

During  the  first  year  of  the  funding  period  we  completed  Aim  1  (Tasks  1-2).  Specifically,  we 
determined  the  2.5A  resolution  structure  of  crystals  that  were  prepared  by  mixing  the  p53  core 
domain  with  the  F1-CDB3  peptide.  Unfortunately,  the  structure  did  not  reveal  ordered  electron 
density  for  the  peptide.  Subsequent  experiments  involved  soaking  preformed  p53  core  domain 
crystals  with  peptide  which  also  produced  a  structure  in  which  no  ordered  density  for  the  peptide 
could  be  identified.  We  conclude  that  the  FL-CDB3  peptide  does  not  bind  p53  in  a  unique 
location  and  conformation  and  therefore  that  it  is  not  possible  to  characterize  a  structure  of  a 
p53/FL-CDB3  complex.  This  is  consistent  with  recent  observations  that  have  been  made  by 
Fersht  and  coworkers  (Friedler  et  al.,  2005). 

During  the  second  and  third  years  of  the  funding  period,  we  completed  Aim  2  (Tasks  3-4),  Aim  3 
(Tasks  5-6)  and  Task  7  of  Aim  4.  Specifically,  we  determined  the  1.55  A  crystal  structure  of  the 
mouse  p53  core  domain  with  a  molecule  of  tris(hydroxymethyl)aminomethane  (Tris)  bound 
through  multiple  hydrogen  bonds  to  a  region  of  p53  shown  to  be  important  for  repair  of  a  subset 
of  tumor-derived  p53  stability  mutations.  Consistent  with  the  hypothesis  that  Tris  binding 
stabilizes  the  p53  core  domain,  we  carried  out  equilibrium  denaturation  experiments  that 
demonstrate  that  Tris  binding  increases  the  thermodynamic  stability  of  the  mouse  p53  core 
domain  by  0.74  kcal/mol;  and  molecular  dynamic  simulations  revealing  an  overall  reduction  in 
root-mean-square  deviations  of  the  core  domain  of  0.7  A  when  Tris  is  bound.  We  also  showed 
that  these  crystals  of  the  p53  core  domain  are  suitable  for  Multiple  Solvent  Crystal  Structure 
approach  to  identify  other  potential  binding  sites  for  possible  core  domain  stabilization 
compounds.  Analysis  of  the  residue-specific  temperature  factors  of  the  high  resolution  core 
domain  structure,  coupled  with  a  comparison  with  other  core  domain  structures,  also  reveals  that 
the  LI,  H1-S5  and  S7-S8  core  domain  loops,  also  shown  to  mediate  various  p53  activities, 
harbor  inherent  flexibility  suggesting  that  these  regions  might  be  targets  for  other  p53  stabilizing 
compounds.  Together,  these  studies  provide  a  molecular  scaffold  for  the  structure-based  design 
of  p53  stabilization  compounds  for  development  as  possible  therapeutic  agents.  These  studies 
are  described  in  our  recent  publication  in  Acta  Crystallographica  (Ho  et  al.,  2006). 

During  the  third  year  of  the  funding  period  we  have  extended  Aim3  (Tasks  5-6).  In  particular, 
we  have  extended  our  findings  that  Tris  stabilizes  the  p53  core  domain  to  identify  Tris-like 
molecules  that  are  predicted  to  bind  p53  with  higher  affinity  and  potentially  to  increase  the 
degree  of  small  molecule  stabilization.  Using  Tris  as  a  lead  compound  for  p53  core  domain 
stabilization,  we  carried  out  a  virtual  screening  using  the  SPECS  (http://www.specs.net/)  and 
TimTec  (http://blaster.docking.org/zinc/)  databases.  Screening  was  performed  on  the  Pittsburg 
Supercomputing  Center  (www.psc.edu)  using  a  Linux  server  in  our  lab.  Since  a  scoring  function 
has  not  yet  been  developed  to  reliably  and  consistently  rank  and  quantitate  ligand-protein 
energies,  a  heuristic  docking  and  consensus  scoring  strategy  was  used  in  the  virtual  screening.  In 
this  particular  case,  the  program  DOCK4.0  (Morris  et  al.,  1998)  was  employed  for  the  primary 
screening  with  a  radius  of  6  A  around  the  Tris  molecule.  During  the  molecular  docking 
calculations,  Kollman-all-atom  charges  were  assigned  to  the  protein,  and  Geisterger-Huckel 
charges  were  assigned  to  tris  molecules  due  to  lack  of  proper  Kollman  charges. 

The  conformational  flexibility  of  the  compounds  from  the  databases  were  considered  in  the 
docking  procedure  and  the  DOCK  suite  was  used  to  evaluate  the  results  using  a  shape  scoring 
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function  and/or  a  function  approximating  the  ligand-receptor  binding  energy.  Following  the 
initial  orientation  and  scoring  evaluation,  a  grid-based  rigid  body  minimization  was  carried  out 
for  the  ligands  to  locate  the  nearest  local  energy  minimum  within  the  receptor  binding  site.  The 
position  and  conformation  of  each  docked  molecule  was  optimized  using  single  anchor  search 
and  a  torsion  minimization  method  in  DOCK4.0.  Fifty  configurations  per  ligand  building  a  cycle 
and  50  maximum  anchor  orientations  were  used  in  the  anchor-first  docking  algorithm.  All 
docked  configurations  were  energy  minimized  using  100  maximum  iterations  and  1 
minimization  cycle. 


Following  molecule  selection  based  on  the  docking  results,  the  top  40000  molecules  from  each 
database  were  selected  for  further  analyses.  These  molecules  were  re-scored  using  the  program 
SLIDE,  XSCORE  and  the  scoring  function  of  AutoDock3.0.  Based  on  the  second  scoring 
results,  13  compounds  (10  from  the  SPECS  database  and  3  from  the  TimTec  database)  were 

selected  for  further  analysis  using 
solution  studies.  Some  representative 
compounds  from  this  set  are  shown  in 
Figure  1.  These  Tris  analogues  are 
currently  being  analyzed  for  p53 
stabilization. 


Figure  1.  Structure  of  4  Tris  like  compounds  that  showed 
favorable  p53  docking  properties. 


In  a  related  study  we  collaborated  with 
the  El-Diery  laboratory  at  the 
University  of  Pennsylvania  to  analyze 
the  properties  of  acridine  molecule,  CP- 
31398,  and  related  compounds  in  p53 
stabilization  and  activation. 

Specifically,  we  showed  that  CP-31398 
activates  wild-type  p53  by  a  novel 
mechanism  that  does  not  involve  phosphorylation  of  the  amino-terminus  of  p53  and 
disassociation  of  MDM2.  To  identify  more  potent  CP-3 1398-like  p53  activators,  we  synthesized 
4  acridine  derivatives  with  a  similar  structure  to  CP-3 1398.  These  four  compounds  induced 
strong  p53  transcription  in  cells  with  wild-type  p53.  We  also  found  that  several  randomly 
chosen  acridine  derivatives,  including  9-  aminoacridine,  amsacrine,  quinacrine  and  acridine 
orange,  induced  p53  transcriptional  activity.  All  these  acridine  derivatives  stabilized  p53  protein 
by  blocking  its  ubiquitination,  without  phosphorylation  of  serl5  or  ser20  on  p53.  Furthermore, 
acridine  derivatives  induced  p53-dependent  cell  death.  Knockout  of  Bax,  a  p53  target  and  a  key 
cell  death  inducer  in  both  intrinsic  and  extrinsic  apoptotic  pathways,  blocked  acridine  derivatives 
from  inducing  cell  death.  In  addition,  in  vivo  delivery  of  quinacrine  and  amsacrine  induced  p53 
transcriptional  activity  in  tumor  xenografts.  Our  results  reveal  that  DNA-intercalating  acridine 
derivatives  can  induce  p53  stabilization  by  a  manner  similar  to  CP-3 1398.  These  findings 
provide  insights  into  p53  regulation  in  response  to  DNA  intercalating  drugs  and  may  assist  new 
anticancer  drug  design.  These  studies  are  described  in  our  recent  publication  in  Cancer  Biology 
and  Therapy  (Wang  et  al.,  2005). 
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(7)  KEY  RESEARCH  ACCOMPLISHMENTS 

We  have  demonstrated,  both  in  solution  and  in  silico,  that  Tris  binding  stabilizes  the  core 
domain  of  p53  and  therefore  Tris  qualifies  as  a  suitable  lead  compound  for  the  structure-based 
optimization  of  p53  stabilizing  compounds  with  possible  therapeutic  application  for  p53-mediate 
breast  cancer. 

We  have  identified  second  generation  Tris-like  p53  stabilizing  compounds  in  silico  that 
are  suitable  for  further  investigation  of  their  p53  stabilization  properties  in  solution. 

We  have  characterized  the  biochemical  properties  of  the  acridine  molecule  CP-3 1398  and 
related  compounds  in  p53  activation  and  stabilization. 

(8)  REPORTABLE  OUTCOMES 

Wang,  W.,  Ho,  W.C.,  Dicker,  D.T.,  Mackinnon,  C.,  Winkler,  J.D.,  Marmorstein,  R.  and  El- 
Deiry,  W.C.  “Acridine  derivatives  activte  p53  and  induce  tumor  cell  death  through  Bax”  (2005) 
Cancer  Biol.  Ther.  44,  893-898. 

Ho,  W.C.,  Luo,  C.,  Zhao,  K.,  Chai,  X.,  Fitzgerald,  M.X.  and  Marmorstein,  R.  “High  resolution 
structure  of  the  p53  core  domain:  implications  for  binding  small  molecule  stabilizing 
compounds.”  (2006)  Acta  Cryst.,  62,  1484-1493. 

(9)  CONCLUSIONS 

In  future  studies,  we  will  carry  out  Aim  4  (Tasks  7-9)  to  further  characterize  the  Tris  analogues 
that  we  identified  through  our  virtual  screening  procedure,  in  solution  for  p53  binding  and 
stability  properties.  Compounds  that  show  increased  solution  binding  and  stability  relative  to 
Tris  will  be  cocrystallized  with  the  p53  core  domain.  We  will  also  continue  our  virtual  screening 
to  identify  additional  Tris  analogues  (hits)  that  are  predicted  to  bind  and  stabilize  the  p53  core 
domain  and  we  will  further  filter  these  hits  for  compounds  that  show  good  solubility 
characteristics.  Additional  "hits”  will  be  further  analyzed  in  solution  as  described  above. 

The  structure-based  drug  design  approach  (often  called  “rational  drug  design”),  that  we  are  using 
towards  the  development  of  small  molecule  compounds  that  might  restore  function  to  tumor- 
derived  p53  mutants,  is  a  recently  exploited  and  particularly  powerful  strategy  which  uses 
protein  structural  information  to  specifically  design  small  peptides  or  non-peptidic  molecules 
that  modulate  the  activity  of  a  protein  of  interest  (Garrett  and  Workman,  1999;  Huang,  2000; 
Jackson,  1997;  Oakley  and  Wilce,  2000;  Tada  et  al.,  1999;  Wang  et  ah,  1999;  Wieczorek  et  ah, 
1996).  This  strategy  has  shown  considerable  promise,  already  yielding  clinically  useful  peptides 
and  compounds  (Amzel,  1998;  Gane  and  Dean,  2000;  Kirkpatrick  et  al.,  1999;  Klebe,  1998; 
Kubinyi,  1998;  Lunney,  1998;  Roe  et  al.,  1998;  Sehgal,  2002)  as  well  as  several  other 
compounds  currently  in  clinical  trials  (Klebe,  1998).  Based  on  our  encouraging  results  to  date, 
we  propose  that  a  structure-based  approach  is  an  effective  strategy  of  achieving  our  ultimate  goal 
of  developing  p53-targeting  drugs  that  will  have  clinical  application  for  the  treatment  of  p53- 
mediated  breast  cancer. 
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High-resolution  structure  of  the  p53  core  domain: 
implications  for  binding  small-molecule  stabilizing 
compounds 


The  p53  transcriptional  regulator  is  the  most  frequently 
mutated  protein  in  human  cancers  and  the  majority  of  tumor- 
derived  p53  mutations  map  to  the  central  DNA-binding  core 
domain,  with  a  subset  of  these  mutations  resulting  in  reduced 
p53  stability.  Here,  the  1.55  A  crystal  structure  of  the  mouse 
p53  core  domain  with  a  molecule  of  tris(hydroxymethyl)amino- 
methane  (Tris)  bound  through  multiple  hydrogen  bonds  to  a 
region  of  p53  shown  to  be  important  for  repair  of  a  subset  of 
tumor-derived  p53-stability  mutations  is  reported.  Consistent 
with  the  hypothesis  that  Tris  binding  stabilizes  the  p53  core 
domain,  equilibrium  denaturation  experiments  are  presented 
that  demonstrate  that  Tris  binding  increases  the  thermo¬ 
dynamic  stability  of  the  mouse  p53  core  domain  by 
3.1  kJ  mor1  and  molecular-dynamic  simulations  are  pre¬ 
sented  revealing  an  overall  reduction  in  root-mean-square 
deviations  of  the  core  domain  of  0.7  A  when  Tris  is  bound.  It  is 
also  shown  that  these  crystals  of  the  p53  core  domain  are 
suitable  for  the  multiple-solvent  crystal  structure  approach  to 
identify  other  potential  binding  sites  for  possible  core-domain 
stabilization  compounds.  Analysis  of  the  residue-specific 
temperature  factors  of  the  high-resolution  core-domain 
structure,  coupled  with  a  comparison  with  other  core-domain 
structures,  also  reveals  that  the  LI,  H1-S5  and  S7-S8  core¬ 
domain  loops,  also  shown  to  mediate  various  p53  activities, 
harbor  inherent  flexibility,  suggesting  that  these  regions  might 
be  targets  for  other  p53-stabilizing  compounds.  Together, 
these  studies  provide  a  molecular  scaffold  for  the  structure- 
based  design  of  p53-stabilization  compounds  for  development 
as  possible  therapeutic  agents. 


Received  28  July  2006 
Accepted  22  September  2006 


PDB  References:  p53-Tris 
complex,  2ioi;  p53— HEPES 
complex,  ioo;  p53-2- 
propanol  complex,  2iom. 


1.  Introduction 

The  p53  tumor-suppressor  protein  is  crucial  in  maintaining 
genomic  integrity.  In  the  event  of  DNA  damage,  p53  regulates 
the  transcription  of  genes  that  lead  to  apoptosis  or  cell-cycle 
arrest  (Polyak  et  a  I. .  1997;  Kastan  et  al. ,  1991).  The  p53  protein 
has  four  domains:  an  amino-terminal  transactivation  domain 
(residues  1^14),  a  core  DNA-binding  domain  (102-292),  a 
tetramerization  domain  (residues  320-356)  and  a  carboxyl- 
terminal  regulatory  domain  (residues  320-356).  Greater  than 
50%  of  human  cancers  are  associated  with  mutations  to  p53 
and  an  estimated  95%  of  all  tumorogenic  mutations  map  to 
the  region  of  the  p53  gene  that  encodes  the  core  domain 
(Bullock  &  Fersht,  2001;  Levine,  1997).  The  missense  muta¬ 
tions  that  occur  in  this  region  have  been  classified  into  two 
groups:  mutations  that  directly  disrupt  p53-DNA  contacts  and 
mutations  that  decrease  the  thermal  stability  of  the  core 
domain. 

Owing  to  the  crucial  role  that  the  core  domain  plays  in  the 
function  of  p53,  it  presents  a  potentially  valuable  therapeutic 
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Table  1 

Data-collection  and  refinement  statistics. 


High  resolution 

2-Propanol  soak 

HEPES 

Space  group 

C2 

C2 

C2 

Unit-cell  parameters 
a  (A) 

92.34 

92.30 

91.57 

6(A) 

44.53 

44.54 

44.71 

c  (A) 

63.02 

62.91 

62.61 

P  (°) 

126.25 

126.23 

125.62 

Resolution  (A) 

30-1.55 

50-2.0 

50-2.0 

No.  of  observations 

154271 

57269 

63041 

No.  of  unique  reflections 

30053 

13740 

13086 

Resolution  range  (A) 

30.0-1.55 

50-2.0 

50-2.02 

Completeness  (%) 

97.4  (99.1) 

96.2  (87.4) 

96.1  (91.9) 

Ila(I) 

31.5  (11.7) 

21.1  (3.92) 

25.8  (5.82) 

^merget 

0.040  (0.128) 

0.063  (0.263) 

0.044  (0.198) 

No.  of  protein  atoms 

1470 

1470 

1470 

No.  of  water  atoms 

256 

159 

176 

No.  of  Zn  atoms 

1 

1 

1 

No.  of  Tris  atoms 

8 

8 

0 

R  factor^  (%) 

18.37 

20.89 

21.00 

*free§  (%) 

22.78 

24.65 

22.47 

Mean  B  value  (A2) 

20.44 

29.43 

33.09 

R.m.s.d.  bonds  (A) 

0.010 

0.011 

0.006 

R.m.s.d.  angles  (°) 

2.20 

1.44 

1.33 

Dihedral  angles  (°) 

27.3 

26.2 

25.4 

Improper  angles  (°) 

1.37 

0.84 

0.88 

t  Emerge  =  E*  E;  14,;  -  <4)1  /Ea  E,  4j-  where  <4>  is  the  mean  intensity  of  symmetry- 
related  reflections.  R  factor  =  Ell^ol  —  4 1 1 / E  141,  where  and  l\  are  the 
observed  and  calculated  structure-factor  amplitudes,  respectively.  §  i?free  is  calculated 
for  5%  of  the  data  that  was  withheld  from  refinement. 

target.  Indeed,  a  significant  number  of  studies  have  addressed 
the  feasibility  of  restoring  function  in  tumor-derived  p53  core¬ 
domain  mutants.  In  the  case  of  DNA-contact  mutants,  intro¬ 
duction  of  amino-acid  residues  that  make  additional  contacts 
to  DNA  have  been  shown  to  rescue  the  function  of  these 
mutants  (Wieczorek  et  al.,  1996).  In  the  case  of  mutations  that 
reduce  the  thermostability  of  the  core  domain,  second-site 
mutations  have  been  identified  which  stabilize  common 
tumor-derived  p53  core-domain  mutants  (Brachmann  et  al ., 
1998).  These  second-site  mutations  have  been  shown  to 
mediate  thermodynamically  favorable  interactions  that 
compensate  for  deleterious  mutations  located  in  other  regions 
of  the  core  domain  (Nikolova  et  al. ,  2000).  Together,  these 
studies  suggest  that  the  development  of  small-molecule 
compounds  that  compensate  for  deleterious  DNA-contact  or 
protein-stability  mutants  may  provide  useful  strategies  to 
restore  function  to  common  tumor-derived  p53  mutant 
proteins. 

Structure-based  drug  design  is  a  powerful  strategy  for  the 
development  of  small-molecule  compounds  to  modulate  the 
function  of  a  protein.  Towards  this  goal,  several  p53  core- 
domain  structures  that  might  be  used  for  structure-based 
design  have  been  reported,  albeit  to  medium  resolution.  The 
X-ray  crystal  structure  of  the  human  core  domain  in  complex 
with  DNA  (Cho  et  al. ,  1994)  and  the  mouse  core  domain 
without  DNA  (Zhao  et  al. ,  2001)  have  been  reported  to 
resolutions  of  2.2  and  2.7  A,  respectively.  In  addition,  the 
structures  of  several  superstable  human  p53  core-domain 
mutants  have  been  determined  (Joerger  et  al .,  2004,  2005)  with 
resolution  ranging  from  1.98  to  1.80  A.  Although  these 


structures  provided  important  insights  into  several  aspects  of 
p53  function,  the  moderate  resolution  of  these  structures 
limits  their  usefulness  in  rational  drug-design  strategies. 

In  an  attempt  to  prepare  a  very  high  resolution  p53  core¬ 
domain  structure  that  is  more  amenable  to  structure-based 
drug  design,  we  crystallized  the  mouse  p53  core  domain  and 
identified  a  new  crystallization  condition  that  produces  crys¬ 
tals  that  diffract  to  1.55  A  resolution  and  we  report  its  refined 
structure  here.  Unexpectedly,  the  structure  reveals  a  bound 
buffer  molecule  of  tris(hydroxymethyl)aminomethane  (Tris) 
and  interestingly  this  Tris  molecule  mediates  multiple 
hydrogen  bonds  to  a  region  of  p53  shown  to  be  important  for 
the  repair  of  a  subset  of  tumor-derived  p53  mutations.  Based 
on  this  observation,  we  performed  equilibrium  denaturation 
experiments  and  molecular-dynamics  simulations  in  the 
presence  and  absence  of  Tris  and  the  results  are  consistent 
with  increased  thermodynamic  stability  of  the  p53  core 
domain  in  the  presence  of  Tris.  We  also  show  that  these 
crystals  of  the  p53  core  domain  are  suitable  for  the  multiple- 
solvent  crystal  structure  approach  to  identify  other  potential 
binding  sites  for  possible  core-domain  stabilization.  Finally,  we 
carried  out  an  analysis  of  local  regions  of  disorder  within  the 
mouse  p53  core  domain  to  provide  insights  into  areas  of 
conformational  flexibility  that  might  be  appropriate  sites  for 
the  binding  of  other  p53-stabilizing  compounds.  Together, 
these  studies  provide  a  molecular  scaffold  for  the  structure- 
based  design  of  p53-stabilization  compounds  for  development 
as  possible  therapeutic  agents. 

2.  Experimental  procedures 

2.1.  Protein  purification  and  crystallization 

The  mouse  p53  core  domain  was  purified  as  described 
elsewhere  (Zhao  et  al .,  2001).  Briefly,  the  pRSETA  (Invi- 
trogen)  bacterial  expression  vector  encoding  the  p53  mouse 
core  domain  (residues  92-292)  was  overexpressed  in  Escher¬ 
ichia  coli  BL21(DE3)  induced  with  0.5  m M  isopropyl  1-thio- 
/S-D-galactopyranoside  (IPTG)  at  a  growth  temperature  of 
288  K.  The  protein  was  purified  using  a  combination  of  cation- 
exchange  (SP-Sepharose)  and  gel-filtration  (Superdex-75) 
chromatography  in  a  buffer  containing  20  nr M  sodium  citrate 
pH  6.1,  100  m M  NaCl,  10  \iM  zinc  acetate  and  1  nr M  TCEP. 
The  protein  was  concentrated  to  ~10  mg  ml  1  using  filtration 
by  centrifugation  and  frozen  at  193  K  for  storage  prior  to 
crystallization.  Crystals  were  obtained  by  the  hanging-drop 
method  by  mixing  protein  and  reservoir  solutions  in  a  1:1  ratio 
and  equilibrating  against  a  reservoir  solution  containing 
100  m M  Tris  pH  7.0  and  16-18%  PEG  2K  MME  (mono- 
methylether)  at  293  K.  Crystals  appeared  overnight  and  grew 
to  typical  dimensions  of  0.3  x  0.1  x  0.1  mm  in  2-3  d.  Crystals 
were  flash-frozen  in  100  m M  Tris  pH  7.0,  20%  PEG  2K  MME 
and  20%  MPD. 

2.2.  Data  collection,  structure  determination  and  refinement 

Two  data  sets  were  collected  at  120  K  at  the  19BM  beamline 
at  the  Advanced  Photon  Source,  Argonne  National  Labora- 
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tory  using  a  Quantum  CCD  detector.  A  high-resolution  data 
set  was  collected  at  a  crystal-to-detector  distance  of  120  mm. 
In  order  to  avoid  overloads  in  the  low-resolution  bins,  a  low- 
resolution  data  set  was  collected  at  a  crystal-to-detector 
distance  of  240  mm  with  a  shorter  exposure  time  using  the 
same  crystal.  Crystals  form  in  space  group  C2,  with  one 
molecule  per  asymmetric  unit,  and  diffract  strongly  to  1.55  A 
resolution  (Table  1).  Data  were  processed  with  HKL- 2000  and 
the  structure  was  solved  by  molecular  replacement  with  the 
program  AMoRe  (Navaza,  1994)  using  data  from  27  to  3.0  A 
resolution  with  the  lhu8  structure  as  a  search  model.  Rotation 
and  translation  searches  followed  by  rigid-body  refinement 
yielded  an  unambiguous  solution  with  an  R  factor  of  35.9% 
and  a  correlation-coefficient  value  of  69.4%.  The  protein 
model  was  initially  refined  automatically  using  data  from  7  to 
1.55  A  resolution  with  the  programs  ARP 1  w ARP  and  guiSIDE 
(Perrakis  et  al. ,  1999;  Perrakis  et  al. ,  2001)  and  manual 
adjustments  were  carried  out  with  the  program  O  (Jones  et  al ., 
1991).  Simulated-annealing  molecular  dynamics  was 
performed  using  CNS  (Briinger  et  al.,  1998).  As  the  refinement 
of  atomic  position  reached  convergence,  isotropic  B  factors 
were  refined  and  water  molecules  were  added  using  the 
waterpick  routine  of  CNS.  The  data  was  then  extended  to  30- 
1.55  A  resolution  and  a  bulk-solvent  correction  was  added. 
SHELX91  was  used  for  further  refinement  (Sheldrick  & 
Schneider,  1997)  against  F  and  the  same  Rbee  set  was  used 
throughout  the  refinement.  The  model  was  subjected  to  a 
stepwise  increase  in  resolution  (STIR)  from  3.0  to  1.55  A  using 
0.005  A  increments,  resulting  in  RWorkmg  and  Rtree  of  21.14  and 
26.01%,  respectively.  Anisotropic  displacement  parameter 
refinement  was  then  introduced  (ANIS)  and  the  model  was 
refined  for  another  30  rounds.  This  resulted  in  a  dramatic  drop 
in  the  Rworidng  and  R(ree  parameters  to  18.83  and  22.90%, 
respectively.  The  resulting  electron-density  maps  were 
considerably  improved.  At  this  point,  clear  electron  density 
could  be  observed  for  a  molecule  of  tris(hydroxymethyl)- 
aminomethane  (Tris),  which  was  in  the  crystallization  buffer. 
Tris  was  then  appropriately  built  using  the  program  O  and 
refined  using  SHELX91,  resulting  in  a  final  model  with 
^working  and  Rtree  parameters  of  18.37  and  22.78%,  respec¬ 
tively.  The  final  model  has  excellent  stereochemistry  (Table  1), 
with  89.1%  (139  total)  of  the  residues  in  the  most  favorable 
regions  and  10.9%  (17  total)  in  additional  allowed  regions  of 
the  Ramachandran  plot. 

p53  core-domain  crystals  for  2-propanol-soaking  experi¬ 
ments  were  harvested  in  100  m M  Tris  pH  7.0,  20%  PEG  2K 
MME  and  soaked  in  100  m M  Tris  pH  7.0,  20%  PEG  2K  MME, 
35%  2-propanol  for  at  least  6  h.  Crystals  were  looped  and 
frozen  directly  in  liquid  propane.  Data  collection  was  carried 
out  on  a  Rigaku  R-AXIS  IV  detector  with  a  rotating  copper- 
anode  generator.  Crystals  diffracted  to  a  maximum  of  2.0  A. 

Eo,  2-propanol  C,  native?  C,  native  E,  2-propanol?  Ei  2-propanol 

—  Fc  2-propanoi  and  2 Fa  —  Fc  maps  were  generated  using  the 
previously  solved  native  structure  with  all  waters  removed. 
Assignment  of  potential  2-propanol  density  was  based  on 
careful  examination  of  electron  density  at  various  a  levels  in 
conjunction  with  plausible  functional  group  interactions.  After 
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addition  of  waters,  the  structure  was  refined  to  final  R(actor  and 
Rfree  parameters  of  20.89  and  24.65%,  respectively. 

Mouse  p53  core-domain  crystals  grown  in  the  absence  of 
Tris  were  obtained  from  conditions  similar  to  the  Tris-bound 
crystals  except  for  the  substitution  of  100  m M  HEPES  pH  7.0 
for  100  nr M  Tris  pH  7.0.  Crystals  appeared  to  be  isonrorphous 
to  those  obtained  using  Tris  buffer  and  were  frozen  using  the 
same  cryoprotectant  conditions.  Data  were  collected  at  the 
Advanced  Photon  Source,  Argonne  National  Laboratory. 
Crystals  diffracted  to  a  maximum  of  2.02  A.  The  unit-cell 
parameters  were  nearly  identical  to  those  of  the  crystals 
obtained  with  Tris  buffer  and  the  model  was  built  using  the 
previously  solved  structure  with  waters  removed  as  a  starting 
structure.  All  structural  figures  were  rendered  using  PyMOL 
unless  otherwise  noted  (DeLano,  2002). 

2.3.  Equilibrium  denaturation 

Protein  was  diluted  from  concentrated  aliquots  to  10  pM 
using  premade  stocks  of  buffer  consisting  of  either  buffers 
A,  B,  C  or  D.  Buffer  A  contained  50  m M  phosphate  pH  6.4, 
100  m M  NaCl  and  1  m M  TCEP  Buffer  B  contained  50  m M 
phosphate  pH  6.4,  100  m M  NaCl,  1  m M  TCEP  in  addition  to 
8  M  urea.  Buffers  C  and  D  were  identical  to  buffers  A  and  B , 
respectively,  except  for  the  presence  of  1  m M  Tris-HCl. 
Measurements  were  taken  on  a  PTI  spectrofluorimeter 
(Photon  Technologies  International,  Birmingham,  NJ,  USA) 
equipped  with  an  LPS-220B  lamp  power  supply,  MD-5020 
motor  driver  and  SC-500  shutter  control  and  operated  using 
the  Felix  software  package.  The  temperature  was  maintained 
at  283  K  using  a  water  bath.  Typically,  800  pi  of  10  p M  protein 
dissolved  in  buffer  A  or  C  was  placed  in  a  quartz  cuvette  with  a 
5  mm  path  length  (Starna  Cells,  Atascadero,  CA,  USA).  An 
excitation  wavelength  of  280  nm  (4  nm  bandpass)  was  used 
and  a  spectrum  was  then  recorded  from  300  to  400  nm  (4  nm 
bandpass).  Fully  unfolded  protein  (8  M  urea)  was  observed  to 
have  a  fluorescence  maximum  at  350  nm.  After  each  reading, 
20  pi  of  the  protein  solution  in  the  cuvette  was  removed  and 
replaced  with  20  pi  of  protein  dissolved  in  either  buffer  B  or  D 
and  then  stirred  using  a  magnetic  stirrer.  The  pH  (6.4)  of  urea 
solutions  at  I/50%  in  the  presence  or  absence  of  Tris  was 
identical.  Unfolding  of  the  mouse  p53  core  domain  in  urea  was 
reversible  as  evidenced  by  the  observation  of  a  similar  spec¬ 
trum  to  folded  protein  upon  dilution  of  the  protein  from  8  M 
urea.  After  normalization,  the  data  were  fitted  to  a  two-state 
equation  (Bullock  et  al.,  1997), 

P  _  (o'a?  +  PN[P])  +  (aD  +  PD[P])  exp {m([D\  -  [D]50%)/RT} 

1  +  exp{m([D]  -  [D]50%)/RT} 

using  the  Origin  software  package  (Microcal,  Northhampton, 
MA,  USA),  where  F  is  the  fluorescence  at  350  nm,  aN  and  aD 
are  the  intercepts  at  fully  folded  and  unfolded  conditions, 
respectively,  fiN  and  jiD  are  the  slopes  at  fully  folded  and 
unfolded  conditions,  [T>]  is  the  concentration  of  urea  and 
[T>]50%  is  the  concentration  of  urea  at  the  point  of  inflection, 
with  m  being  the  slope  of  that  transition.  A  AG  for  Tris  and 
protein  or  protein  alone  were  calculated  using  the  formula 
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A  AG  -  (ttl)  (  D|ns  50%  -^protein  alone,  50%)  (Bullock  et  Cll. ,  1997) 
with  ( m )  values  calculated  as  the  average  m  values  of  the  two 
data  sets  (—2.41  ±  0.12  and  —2.20  ±  0.08,  respectively). 

2.4.  Molecular-dynamics  simulation 

Initial  coordinates  for  the  p53-Tris  complex  and  p53  (Tris 
removed  from  the  p53-Tris  complex)  were  taken  from  the 
high-resolution  crystal  structures.  All  molecular-dynamics 
(MD)  simulations  were  performed  using  the  parallelized  MD 
program  GROMACS  3.3  with  constant  number,  pressure  and 
temperature  (NPT)  and  periodic  boundary  conditions  (Stryer, 
1988;  Kabsch  &  Sander,  1983;  Mierke  &  Kessler,  1991).  Before 
MD  simulations,  explicit  polar  and  aromatic  H  atoms  were 
added  for  both  models.  The  GROMOS96  force  field  was  used 
for  the  proteins.  The  topology  file  and  other  force-field  para¬ 
meters  of  Tris  were  generated  using  the  PRODRG  program 
(Schiittelkopf  &  van  Aalten,  2004). 


Both  the  p53-Tris  complex  and  p53  were  solvated  in  a 
rectangular  periodic  box  filled  with  SPC  water  molecules 
(Berendsen  et  al.,  1981).  The  minimum  distance  between  the 
protein  atoms  and  the  box  walls  was  set  to  be  >10  A. 
Counterions  were  added  to  neutralize  the  net  charges  of  the 
systems  and  each  system  was  subjected  to  energy  minimization 
using  the  steepest-descent  method  until  the  energy  conver¬ 
gence  at  100  kJ  moP1  nnP  .  Afterwards,  each  system  was 
equilibrated  for  200  ps.  Finally,  a  5  ns  MD  simulation  was 
performed  on  each  system.  The  linear  constraint  solver 
(LINCS)  method  (Hess  et  al.,  1997)  was  used  to  constrain 
bond  lengths,  allowing  an  integration  step  of  2  fs.  Electrostatic 
interactions  were  calculated  with  the  Particle-Mesh  Ewald 
algorithm  (Darden  et  al.,  1993;  Essmann  et  al.,  1995).  A 
constant  pressure  of  10s  Pa  was  applied  with  a  coupling 
constant  of  1.0  ps.  The  Particle-Mesh  Ewald  (PME)  method 
was  used  to  calculate  electrostatic  interactions.  The  simulated 
systems  were  coupled  into  an  external  temperature  bath  at 
300  K  with  a  coupling  constant  of  0.1  ps  and  isotropic  pressure 
coupling  with  time  constant  of  1  ps  was  applied  to  keep  the 
pressure  at  10s  Pa.  The  structures  for  analysis  were  saved 
every  500  steps  (1  ps).  MD  simulations  were  run  on  Pittsburgh 
Supercomputing  Center  using  Linux  servers.  Analyses  were 
performed  using  facilities  within  the  GROMACS  package. 

3.  Results 

3.1.  Overall  structure  of  the  mouse  p53  core  domain 

The  crystal  structure  of  the  mouse  p53  core  domain  was 
determined  to  1.55  A  in  space  group  C2  and  represents  the 
highest  resolution  structure  of  the  core  domain  to  date.  The 


(c) 


Figure  1 

Overall  structure  of  the  mouse  p53  core  domain,  (a)  Representative  F0  —  Fc  OMIT  electron-density  map  of  residues  231-233  contoured  at  3.5cr.  ( b ) 
Overall  structure  of  the  mouse  p53  core  domain  solved  to  1.55  A  resolution.  The  Zn  atom  is  shown  as  a  yellow  sphere  and  a  molecule  of  Tris  is  shown  in 
red.  (c)  Overlay  of  the  mouse  core  domain  solved  at  1.55  A  resolution  in  space  group  Cl  (HR,  green),  the  mouse  p53  core  domain  solved  at  2.7  A  in 
space  group  C2221  (PDB  code  lhu8,  grey)  and  the  human  p53  core  domain  unbound  to  DNA  (PDB  code  ltsr;  orange) 
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p53  core  domain  crystallizes  as  a  monomer  and  contains  one 
molecule  per  asymmetric  unit.  The  electron-density  map  is  of 
excellent  quality  (Fig.  la),  with  virtually  all  of  the  main  chain 
and  a  majority  of  the  side  chains,  including  surface  lysine  and 
methionine  residues,  which  are  typically  disordered,  modeled 
unambiguously  using  2 Fa  —  Fc  and  Fa  —  Fc  electron  density. 


The  overall  fold  of  the  high-resolution  mouse  p53  core 
domain  (Fig.  lb)  is  similar  to  that  previously  reported  for 
other  p53  core-domain  structures  at  lower  resolution  (Zhao  et 
al.,  2001;  Cho  et  al.,  1994;  PDB  codes  lhu8  and  ltsr,  respec¬ 
tively).  Briefly,  the  core  domain  forms  two  regions  of  twisted 
antiparallel  sheets  (SI,  S3,  S5,  S8  and  S6,  S7,  S4,  S9,  S10).  In 


Residue  No. 
(e) 


Figure  2 

(a)  Fa  —  Fc  omit  density  of  Tris  bound  to  the  p53  core  domain  contoured 
at  2.0cr.  Hydrogen-bonding  interactions  are  depicted  in  orange.  Two 
molecules  of  water  are  depicted  in  red.  ( b )  2 Fa  —  Fc  map  (orange) 
contoured  at  2.5cr  around  the  Tris-binding  region  of  the  mouse  p53  core¬ 
domain  structure  in  the  absence  of  Tris.  Overlaid  as  blue  density  is  an 
Fa  —  Fc  OMIT  map  calculated  at  2.5a  with  W53  omitted.  Crystals  were 
obtained  in  100  mM  HEPES  pH  7.0,  16-18%  PEG  2K  MME 
(monomethylether)  at  293  K.  (c)  Fluorescence  of  the  mouse  p53  core 
domain  at  350  nm  as  a  function  of  urea  concentration.  Normalized  and 
averaged  denaturation  curves  are  shown  in  red  and  back  for  protein  in 
the  absence  and  presence  of  1  m M  Tris,  respectively.  Standard  deviations 
from  the  mean  are  shown  as  error  bars  (n  =  2  for  protein  without  Tris, 
n  =  4  for  protein  with  1  m M  Tris  added),  (rf)  The  r.m.s.d.  values  for  the 
p53-Tris  complex  and  p53  from  molecular-dynamics  simulations,  (e) 
Residue-specific  r.m.s.d.  values  for  the  p53-Tris  complex  and  p53  from 
molecular-dynamics  simulations. 
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addition,  there  are  two  a-helices:  one  at  the  C-terminal  end  of 
the  molecule  (H2)  and  another  shorter  helix  (HI)  near  a  site 
for  zinc  chelation.  A  structural  Zn  atom  is  coordinated  by  four 
residues:  cysteines  238,  242  and  176  and  histidine  179  (human 
numbering).  Also  included  in  the  high-resolution  structure  is  a 
well  ordered  Tris  buffer  molecule  that  is  bound  on  the  protein 
surface  between  the  SI  and  S10  strands  (Fig.  lb). 

3.2.  Comparison  with  other  p53  core-domain  structures 

The  high-resolution  mouse  p53  core-domain  structure  in 
space  group  C2  (HR)  superimposes  well  with  the  2.7  A 
structure  of  the  mouse  p53  core  domain  in  space  group  C222 , 
(PDB  code  lhu8)  and  both  the  DNA-bound  and  unbound 
forms  (chains  B  and  A  of  PDB  code  ltsr)  of  the  2.2  A  human 
p53  core-domain  structure  (Cho  et  al.,  1994).  A  superposition 
of  the  structures  reveals  r.m.s.  deviation  values  of  0.447,  0.753 
and  0.670  A  for  C"  atoms,  respectively.  Not  surprisingly,  the 
central  /1-sheet-rich  core  region  shows  the  greatest  structural 
superposition,  while  the  surrounding  loops  show  the  greatest 
divergence.  Specifically,  the  LI  (between  the  SI  and  S2 
strands)  and  H1-S5  loops  show  the  greatest  structural  diver¬ 
gence.  This  result  highlights  the  high  degree  of  structural 
conservation  between  the  core  domains  of  p53  from  different 
species  and  reveals  that  DNA-specific  binding  does  not 
involve  significant  protein  structural  rearrangement,  which  is 
consistent  with  previous  proposals  (Zhao  et  al. ,  2001). 

A  comparison  of  the  DNA-bound  form  with  HR  reveals 
that  while  the  H2  helix,  jointly  responsible  for  DNA  binding 
with  the  LI  loop,  remains  essentially  unchanged,  the  LI  loop 
undergoes  a  larger  structural  shift  between  the  two  proteins 
(Fig.  lc).  This  result  was  not  surprising  since  the  LI  loop  is 
apparently  flexible  as  revealed  by  comparison  of  several  p53 
core-domain  crystal  structures  (PDB  codes  2bin,  2bio,  2bip, 
2biq,  2geq).  Additionally,  residues  178-186  of  the  mouse  p53 
core  domain,  corresponding  to  the  H1-S5  loop,  shows  a  C“ 
r.m.s.  deviation  of  0.438  A  when  compared  with  the  lower 
resolution  mouse  p53  core  domain  (PDB  code  lhu8).  The 
lhu8  structure  implicates  the  H1-S5  loop  of  one  p53  core 
domain  in  mediating  dimerization  contacts  via  interactions 
with  the  S4-H1  (L2)  and  S6-S7  loops  of  another  p53  core 
domain  (Zhao  et  al.,  2001).  This  interaction  maintains  p53  in  a 
configuration  that  is  incompatible  with  simultaneous  DNA 
binding  by  both  protomers  of  the  dimer.  In  HR,  the  H1-S5 
loop  is  also  involved  in  intersubunit  crystal  contacts,  but  with 
the  S5-S6  loop.  This  comparison  suggests  that  the  flexibility  of 
the  H1-S5  loop  could  play  an  important  role  in  p53  function 
via  intersubunit  core-domain  interactions. 

3.3.  Location  and  characterization  of  a  Tris  molecule  in  the 
core-domain  structure 

A  molecule  of  Tris  was  observed  in  the  electron-density 
map  of  HR  during  refinement  of  the  model.  The  electron 
density  for  this  molecule  is  very  well  defined  and  the  average 
thermal  factor  for  the  Tris  atoms  is  32.50  A2,  which  is  within 
the  average  range  for  thermal  factors  for  the  HR  structure. 
The  Tris  molecule  binds  between  the  SI  and  S10  strands  of  the 


protein  and  makes  direct  and  water-mediated  hydrogen¬ 
bonding  interactions  with  the  side-chain  atoms  of  residues 
Tyrl23,  Asnl28  and  Asp265  (Fig.  2a).  Additionally,  the  Tris 
molecule  makes  water-mediated  hydrogen-bonding  interac¬ 
tions  to  the  backbone  carbonyl  of  Ser266  and  to  the  side-chain 
N  atom  of  GlnlOl.  Interestingly,  the  location  of  the  bound  Tris 
molecule  is  close  to  the  position  of  the  human  N268D  stabi¬ 
lization  mutation  and  appears  to  play  a  similar  role  in  intro¬ 
ducing  bridging  interactions  between  the  SI  and  S10  strands  of 
the  p53  core  domain. 

To  further  assess  the  structural  significance  of  the  Tris- 
mediated  p53  interaction,  we  prepared  isomorphous  crystals 
of  the  mouse  p53  core  domain  in  which  HEPES  buffer  was 
used  in  place  of  Tris  buffer  and  we  determined  the  structure  to 
2.0  A  resolution  (Table  1).  The  structure  of  the  p53  core 
domain  in  the  presence  of  HEPES  buffer  reveals  that  one 
water  molecule  (W53)  overlaps  a  region  of  electron  density 
that  had  been  occupied  by  Tris  and  each  of  the  two  water 
molecules  that  had  bridged  Tris-p53  interactions  (W27  and 
W36)  are  also  present  (Fig.  2b).  Notably,  however,  W53  does 
not  appear  to  compensate  for  the  Tris-mediated  intra-protein 
interactions. 

Although  attempts  to  observe  direct  binding  of  Tris  to  the 
p53  core  domain  in  solution  using  ITC  (isothermal  titration 
calorimetry)  were  unsuccessful  (RM  and  WCH,  unpublished 
data),  presumably  owing  to  the  weak  binding  affinity  of  Tris 
for  the  p53  core  domain,  we  tried  to  test  directly  whether  high 
concentrations  of  Tris  (1  m M)  could  stabilize  the  p53  core 
domain  by  carrying  out  chemical  denaturation  experiments  of 
the  p53  core  domain  in  the  presence  or  absence  of  Tris.  For 
these  studies,  the  mouse  p53  core  domain  was  subjected  to 
urea  equilibrium  denaturation  experiments  in  the  presence  or 
absence  of  1  m M  Tris.  Unfolding  was  monitored  by  measuring 
fluorescence  of  a  single  buried  tryptophan,  which  increased  in 
fluorescence  to  a  maximum  of  X  =  350  nm  as  unfolding 
proceeded.  These  experiments  revealed  that  in  the  presence  of 
1  m M  Tris,  the  mouse  p53  core  domain  has  a  U50%  (concen¬ 
tration  of  urea  required  for  50%  unfolding)  of  2.58  ±  0.02  M, 
compared  with  2.20  ±  0.08  M  for  protein  in  the  absence  of  Tris 
(Fig.  2c).  Analysis  of  this  data  yields  a  AA G  of  3.1  kJ  mol  1 
for  Tris  binding. 

3.4.  Molecular-dynamic  simulations  of  the  p53  core  domain 
in  the  absence  and  presence  of  bound  Tris 

The  behavior  of  the  p53-Tris  complex  and  p53  alone  was 
studied  by  molecular-dynamics  simulation  to  account  for 
protein  flexibility  within  a  simulated  solution  environment. 
The  starting  structures  extracted  from  the  p53  and  p53-Tris 
crystal  structures  reported  here  were  subjected  to  5.0  ns  MD 
simulations  and  the  r.m.s.d.  values  of  backbone  atoms  from 
their  initial  positions  (f  =  0  ps)  were  used  to  measure  core¬ 
domain  stability  and  to  gain  insight  into  possible  structural 
fluctuation.  The  time  evolution  of  the  backbone  r.m.s.d.  values 
of  p53  core-domain  atoms  for  both  systems  is  presented  in 
Fig.  2(d).  In  the  plot,  a  sharp  rise  in  r.m.s.d.  was  observed 
during  the  first  200  ps  for  all  residues  and  was  followed  by  a 
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relative  flattening  out  of  these  fluctuations.  However,  the 
magnitude  of  these  r.m.s.d.  curves  does  not  continue  to 
increase  after  about  a  1.5  ns  MD  simulation,  implying  that 
both  systems  are  stable  over  this  timescale.  The  average 
r.m.s.d.  values  are  about  3.0  and  3.7  A,  respectively,  for  the  p53 
core  domain  in  the  presence  and  absence  of  bound  Tris. 
Interestingly,  the  reduced  r.m.s.d.  values  in  the  presence  of  Tris 
are  not  localized  but  are  distributed  throughout  the  p53  core¬ 
domain  structure  (Fig.  2e). 

3.5.  Multiple-solvent  crystal  structure  analysis  using  p53 
core-domain  crystals 

Although  we  had  fortuitously  discovered  Tris  to  be  a  p53- 
stabilizing  compound,  a  more  systematic  approach  for  deter¬ 
mining  possible  binding  sites  on  the  surface  of  a  protein  which 
can  potentially  be  targeted  by  small  molecules  is  to  soak  high 
concentrations  of  organic  solvents  into  the  crystal  lattice  and 
locate  them  in  high-resolution  electron-density  maps  (English 
et  al.,  1999,  2001;  Mattos  &  Ringe,  1996).  Common  solvents 
used  for  this  purpose  include  2-propanol,  phenol,  dioxane, 
hexane,  acetonitrile  and  acetone,  all  of  which  represent  probes 
for  varying  functional  groups.  The  well  diffracting  HR  crystals 
of  the  mouse  p53  core  domain  were  suitable  for  this  technique. 
Using  this  strategy,  we  soaked  the  HR  crystals  individually 
with  the  compounds  phenol,  2-propanol,  acetonitrile  and 
acetone.  Of  the  molecules  tested,  soaking  p53  crystals  in  35% 
2-propanol  for  6  h  produced  crystals  that  did  not  show 
compromised  diffraction  properties  and  upon  structure 
determination  revealed  the  presence  of  a  2-propanol-binding 
site  within  the  core  domain  (Fig.  3a).  Specifically,  after  data 
collection  and  examination  of  the  electron  density  using  a  com¬ 
bination  of  F0  2-propanol  native?  ^o,  2-propanol  ^*c,  2-propanol 

and  2 Fa  —  Fc  maps,  we  located  one  2-propanol-binding  site 
near  the  LI  loop  of  the  protein  and  bridging  interactions 

between  the  SI  and  S3  /1-sheets.  OMIT  maps  of  the  final 

structure  with  the  2-propanol  omitted  were  generated  to 
confirm  the  presence  of  the  2-propanol  molecule  (Fig.  3b). 

A  high-resolution  view  of  the  interactions  made  between 
2-propanol  and  the  p53  core  domain  reveals  that  the  hydroxyl 
group  of  the  2-propanol  hydrogen  bonds  to  the  backbone 
carbonyl  of  Prol39  and  the  backbone  amide  N  atom  of 

Leulll.  Analysis  of  F0>  2-proPanoi  -  -Fo,  native  as  well  as 

Fa,  native  -  -P0,  2-propanoi  difference  density  at  the  2-propanol- 
binding  site  reveals  several  water  molecules  that  are  displaced 
and  an  altered  hydrogen-bonding  network  in  order  to 
accommodate  the  2-propanol  molecule.  Difference  density 
also  shows  that  the  sulfur  of  Cysl21  also  shifts  conformations 
to  accommodate  the  2-propanol.  The  bottom  of  this  pocket  is 
lined  with  hydrophobic  residues,  including  PhellO,  Phe267, 
Leul30  and  Tyrl23,  which  are  in  a  suitable  environment  to 
interact  with  the  methylene  groups  of  2-propanol.  A  surface 
plot  of  the  region  reveals  a  depression  where  the  2-propanol 
molecule  binds  (Fig.  3c),  which,  interestingly,  is  located  at  a 
one  potential  key  area  of  function  of  the  p53  core  domain,  the 
LI  loop  that  can  participate  in  DNA  binding.  Interestingly,  the 
presence  of  the  2-propanol  compound  at  the  LI  DNA-binding 


loop  suggests  that  this  binding  site  might  provide  a  useful 
template  for  increasing  the  stability  of  the  p53  core  domain 
and  potentially  also  for  introducing  additional  p53  core 
domain-DNA  contacts  that  might  compensate  for  tumor- 
derived  p53  mutations  that  directly  disrupt  p53-DNA 
contacts. 

3.6.  Analysis  of  local  regions  of  disorder  in  the  mouse  p53 
core  domain 

As  the  data-to-parameter  ratio  increases  at  higher  resolu¬ 
tion,  one  can  begin  to  examine  the  regions  of  disorder  of  the 
mouse  p53  core  domain  in  detail,  which  can  provide  further 
insight  into  areas  of  conformational  flexibility  that  might  be 
targets  for  p53-stabilizing  compounds.  To  carry  out  this 
analysis  using  the  high-resolution  mouse  p53  core-domain 
data,  isotropic  thermal  factors  ( B  values)  and  anisotropic 


C b )  (c) 

Figure  3 

Location  of  a  2-propanol  molecule  bound  to  the  mouse  p53  core  domain. 
(a)  Density  surrounding  the  location  of  the  2-propanol-binding  site.  Blue 
density,  2-ProPanoi  -  To,  native  at  the  2.5ff  level,  (b)  Fa  -  Fc  OMIT  map 
calculated  at  2.5 a  using  the  final  structure  with  2-propanol  omitted.  Red 
density,  Fa  native  —  Fa  2_propanob  blue  spheres,  water  molecules  observed  in 
native  structure;  green  sphere,  water  observed  in  2-propanol-soaked 
structure;  orange  dashes,  hydrogen  bonding,  (c)  Surface  representation  of 
the  2-propanol-binding  site  with  2-propanol  shown  as  a  stick  figure.  The 
2-propanol-binding  site  is  denoted  by  an  arrow.  C,  O  and  N  atoms  are 
colored  green,  red  and  blue,  respectively. 
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displacement  parameters  were  refined  as  the  data  allowed 
(Fig.  4a).  A  plot  of  the  average  main-chain  thermal  factors  as  a 
function  of  main-chain  atoms  reveals  that  the  LI,  H1-S5  and 
S7-S8  loops  have  the  highest  thermal  factors,  consistent  with 
the  highest  regions  of  p53  core-domain  flexibility  as  suggested 
from  the  comparison  between  different  p53  core  domains 
(Fig.  lc). 

Residues  within  the  S7-S8  loop  display  the  highest  thermal 
factors  in  the  entire  structure  (Fig.  4b,  top  panel).  Moreover, 
the  main-chain  atoms  of  the  S7-S8  loop  region  also  display  the 
highest  degree  of  anisotropy  in  the  entire  protein  (Fig.  4b, 
bottom  panel).  These  findings  suggest  that  the  S7-S8  loop  is 
not  only  flexible,  but  also  contains  a  directional  component  to 
that  flexibility.  A  plot  of  the  main-chain  thermal  ellipsoids  with 
a  20%  probability  sphere  reveals  a  possible  concerted  ‘hinge- 
like’  movement  of  the  S7-S8  loop  (Fig.  4 d).  On  visual 
inspection,  the  anisotropic  disorder  in  the  main  chain  for 
residues  221,  222,  223  and  224  point  in  the  same  direction, 
possibly  creating  a  hinge  that  pivots  between  Pro220  at  the 


N-terminal  end  and  Glu225  at  the  C-terminal  end.  Analysis  of 
the  Rosenfield  matrix  shows  that  residues  in  the  S7-S8  loop  do 
indeed  move  in  a  direction  independent  from  the  rest  of  the 
protein  (Fig.  4c).  The  high  degree  of  sequence  conservation 
among  p53  proteins  in  the  S7-S8  loop  suggests  that  this  hinge¬ 
like  motion  may  have  a  functional  consequence,  which  may  be 
to  facilitate  intersubunit  contacts.  This  result  is  supported  by 
our  recent  modeling  of  a  tetrameric  complex  from  the  struc¬ 
ture  of  dimeric  mouse  p53  core  domain  bound  to  DNA  (Flo  et 
al,  2006).  In  this  model,  we  predict  that  in  the  tetrameric  form 
the  S7-S8  and  L2  loops  of  the  p53  core  domain  are  likely  to 
make  tail-tail  interactions  when  bound  to  DNA. 

The  LI  loop,  which  can  be  involved  in  DNA  binding,  shows 
the  greatest  degree  of  core-domain  variability  as  determined 
by  the  elevated  thermal  factors  calculated  for  this  region  and 
as  determined  by  comparison  to  the  other  p53  core-domain 
structures  (Fig.  lc).  Interestingly,  analysis  of  the  mean  aniso¬ 
tropy  of  this  region  of  the  structure  reveals  that  the  disorder  is 
more  isotropic.  Further  analysis  of  thermal  ellipsoids  supports 


Figure  4 

Mobility  of  the  p53  core  domain  (a)  Plot  of  temperature  factors  (B)  on  the  secondary  structure  of  the  mouse  p53  core  domain.  Colors  represent 
temperature-factor  magnitudes.  The  highest  temperature  factors  are  in  red  and  the  lowest  are  in  darker  shades  of  blue.  The  figure  was  produced  using 
SWISS-PDB  VIEWER  (Guex  &  Peitsch,  1997)  and  then  rendered  using  POV-Ray.  ( b )  Top,  average  temperature  factor  (B)  of  main-chain  atoms  as  a 
function  of  residue  number;  bottom,  mean  anisotropy  of  main-chain  atoms  plotted  as  a  function  of  residue  number.  Anisotropy  is  shown  as  a  ratio  of  the 
degree  of  displacement  of  the  shortest  axis  to  the  displacement  of  the  longest  axis.  The  figures  were  produced  using  PARVAT1  (Merritt,  1999).  (c) 
Rosenfeld  matrix  generated  using  the  program  ANISOANAL.  Sections  of  similar  shading  indicate  possible  rigid  bodies,  (d)  Thermal  ellipsoid  plot  of 
20%  probability  spheres  for  main-chain  atoms  in  the  S7-S8  loop  region.  Atoms  colored  blue  represent  a  minimal  B  factor  of  10  A2;  atoms  colored  red 
represent  a  maximum  B  factor  of  65  A2.  The  figure  was  produced  using  RASTEP  (Merritt  &  Bacon,  1997).  (e)  Thermal  ellipsoid  plot  of  20%  probability 
sphere  of  the  LI  region.  The  figure  was  produced  using  RASTEP  (Merritt  &  Bacon,  1997). 
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this  (Fig.  4e),  as  many  ellipsoids  are  not  as  elongated  as 
observed  in  the  case  of  the  S7-S8  loop.  The  isotropic  nature  of 
the  disorder  suggests  that  while  the  LI  loop  is  very  flexible,  the 
flexibility  does  not  favor  any  particular  direction.  Analysis  of 
the  main-chain  thermal  factors  of  the  LI  loop  from  the  human 
p53  core  domain-DNA  complex  reveals  thermal  factors  that 
are  about  average  for  the  overall  structure  (Cho  et  al.,  1994), 
reinforcing  the  conclusion  that  the  LI  loop  can  become 
ordered  upon  DNA  binding  by  p53. 

4.  Discussion 

In  this  study,  we  report  a  1.55  A  resolution  structure  of  the 
mouse  p53  core  domain,  representing  the  highest  resolution 
p53  core-domain  structure  reported  to  date.  A  comparison  of 
the  high-resolution  p53  core-domain  structure  reported  here 
with  other  p53  core-domain  structures  from  different  species 
and  in  different  crystal  lattices  reveals  a  highly  conserved  core 
structure  and  three  surface  loops  that  appear  to  be  particularly 
flexible.  These  are  the  S7-S8,  H1-S5  and  LI  loops.  Although 
the  functional  importance  of  the  flexibility  of  the  S7-S8  loop  is 
not  yet  clear,  it  appears  that  the  flexibility  of  the  LI  and  H1-S5 
loops  are  functionally  important.  In  particular,  the  H1-S5  loop 
has  been  proposed  to  play  a  role  in  modulating  the  dimer¬ 
ization  properties  of  the  p53  core  domain  for  suitable  DNA 
binding  and  the  appropriate  conformation  of  the  LI  loop  may 
be  important  for  mediating  sequence-specific  interactions. 

The  high-resolution  p53  core-domain  structure  fortuitously 
uncovered  a  stabilizing  Tris  buffer  molecule  between  the  SI 
and  S10  strands  of  the  core  domain.  Interestingly,  two  Tris 
molecules  were  also  found  in  similar  positions  in  the  recently 
reported  dimeric  mouse  p53-DNA  complex  (Ho  et  al. ,  2006). 
The  Tris  molecule  makes  numerous  interactions  with  the 
protein  that  are  reminiscent  of  the  stabilizing  interactions 
made  by  the  aspartic  acid  of  a  human  N268D  mutant  that 
restored  activity  to  several  tumor-derived  p53  mutant  proteins 
(Joerger  et  al.,  2004).  Correlating  with  the  significance  of  the 
p53-Tris  interactions  seen  in  the  crystals,  equilibrium  dena- 
turation  experiments  demonstrate  that  Tris  increases  the 
thermostability  of  the  p53  core  domain  by  about  3.1  kJ  mol  1 
and  molecular-dynamics  simulations  showed  larger  r.m.s.d. 
values  when  Tris  was  not  bound.  To  our  knowledge,  this 
represents  the  first  identification  of  a  small-molecule  (low- 
molecular-weight)  compound  that  interacts  with  and  increases 
the  thermodynamic  stability  of  p53  and  we  propose  that  the 
p53-Tris  complex  may  provide  a  useful  scaffold  for  the 
structure-based  design  of  p53-stabilizing  compounds. 

The  multiple-solvent  crystal  structure  strategy  of  soaking 
the  crystal  lattice  with  organic  solvents  was  also  employed 
with  the  solvent  2-propanol,  revealing  an  area  of  weak  affinity 
located  at  an  area  of  functional  significance  as  well  as 
considerable  flexibility:  the  LI  loop.  The  2-propanol-binding 
site  might  also  be  exploited  for  the  development  of 
compounds  that  may  stabilize  the  p53  core  domain  and/or 
introduce  additional  p53  core  domain-DNA  interactions  that 
might  compensate  for  a  subset  of  tumor-derived  p53  muta¬ 
tions. 
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Previous  reports  have  described  the  preparation  of  small- 
molecule  compounds  that  restore  p53  function  to  common 
tumor-derived  p53  mutant  proteins.  These  include  a  group  of 
acridine  molecules  (Foster  et  al.,  1999)  and  PRIMA-1  (Bykov 
et  al.,  2002).  However,  more  recent  studies  have  revealed  that 
the  acridine  molecules  do  not  bind  directly  to  p53  (Rippin  et 
al.,  2002).  In  addition,  we  have  recently  found  that  the  PRIMA 
compound  also  does  not  bind  directly  to  p53  (unpublished 
data)  and  thus  is  proposed  to  promote  p53  function  indirectly. 
In  contrast  to  these  earlier  studies,  the  studies  reported  here 
provide  a  rational  framework  for  the  design  and  development 
of  small-molecule  compounds  that  could  directly  target  tumor- 
derived  p53  mutant  proteins  for  the  treatment  of  human 
cancer. 
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